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1. INTRODUCTION 


Life support systems for manned space missions use a combination of replen- 
ishment and recycle. As mission durations increase and crew sizes enlarge, 
either partial or total recycle options become more attractive (Model 1 1977). 

Potential life support systems that provide the essentials- of air, water, 
and food through waste recycling have been referred to as CELSS (Closed Ecologi- 
cal Life Support Systems) and PCELSS (Partially Closed Ecological LSS) 

(Meissner, Model 1 1979). 

For each possible partial or total recycle scenario a break-even point 
exists at which it becomes more economical to send into space all the equipment 
and materials needed to close the recycle loops than to stock all the needed 
provisions for use on a once through basis (Model!, Spurlock 1979). Presently 
the knowledge needed to identify the optimum recycle PCELSS or CELSS option 
for a given type of mission does not exist. 

A scenario analysis method has been proposed for the initial step of 
comparing scenario options and identifying promising alternatives (Modell, 
Spurlock 1979). The method consists of five steps: 

1) Specify the diet. 

2) Select the food-producing processes. 

3) Determine waste-processing requirements and methods. 

4) Characterize equipment and facility requirements. 

5) Test control strategies and determine the ability of the scenario 
to ensure man's survival. 
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This paper presents the results of applying the scenario analysis method 
on a simplified CELSS scenario. Emphasis is on the fifth step of^the method. 
Control strategies and survivability are evaluated with a new approach to the 
analysis of environmental systems developed by Hornberger and Spear (.1980). 

The approach combines probabilistic Mbnte-Carlo simulation techniques with the 
notion of descriptors of system behavior to evaluate system performance. The 
ability of the simplified scenario to ensure man's survival is investigated 
with this approach. The approach can also be used as a generalized sensitivity 
analysis procedure to isolate the critical uncertainties in the system and 
derive information of use in focusing the next phase of research. 

The results presented can not be viewed as a prediction of the survivabil- 
ity of any future CELSS. The scenario model developed is simplified to such 
an extent that it does not have that kind of predictive value. Much of the 
information needed to develop even first preliminary predictive models has not 
yet been collected. 

Instead, the model has been developed and the scenario analysis approach 
applied to it to serve as an example of how the approach may be utilized on 
more realistic predictive models to evaluate the attractiveness of alternative 
partial or total recycle systems. One goal of this paper is to create an 
awareness of the Monte-Carlo System Descriptor evaluation process among those 
responsible for supplying information to the system modelers. 

2. CLOSED ECOSYSTEMS 

A closed ecosystem can be considered as a collection of living and non- 
living components combined in such a way that they are closed to mass exchange ' 
across the total system boundary. Iv'ithin the boundary, however, mass may 
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collect in a flow among the components as determined by the dynamics of the 
ecosystem. 

To survive indefinitely, the mass flows must form closed cycles. The 
simplest ecosystem cycle is a three step process (Quinlan 1975, Quinlan, Paynter 
1975) illustrated in Figure 1. First, plants (autotrophs) use radiant energy 

to synthesize complex organics from carbon dioxide and other nutrients. Next, 

\ 

the plants are eaten by food consumers (.heterotrophs). Ultimately the con- 
sumer's oxidized v^aste products return to the inorganic nutrient pool to com- 
plete the cycle. 

A simple CELSS that includes closed mass flow cycles is shown conceptually 
in Figure 2. In this envisioned space colony humans are the heterotrophs. 

The autotrophs are the plants required to provide the food supply. The human's 
metabolism oxidizes to carbon dioxide and water a portion of the food. The 
remainder appears as partially oxidized products in the human's v;aste products. 
If the recycler fully oxidizes the partially oxidized wastes and any nonedible 
portion of the plants, the inorganic nutrients needed for plant photosynthesis 
are reproduced. The net effect of human metabolism plus v/aste oxidation is 
the reverse of the net plant photosynthesis and the gross material balance 
cycle is closed (Model!, Spurlock 1979). 

Note that the CELSS shown in Figure 2 contains two basic types of elements 
— flow elements and storage elements. The flow elements are considered to 
be "transducers" that transduce a stream of inputs to a steady stream of out- 
puts. With the exception of the growing food, a negligible fraction of the 
system mass is contained in the flow elements. Most mass is held in the vari- 
ous storages. 
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3. MODEL DEVELOPMENT ' ' 

A state equation model based on the conceptual CELSS of Figure 2 has been 
developed to demonstrate scenario and control strategy evaluation by means of 
Monte-Carlo simulation v/ith descriptors of system survival. The five, step 
scenario analysis method guided the model development. 

The inputs required to support a person in space are shown in Table 1. 

A simple diet scenario based on wheat can be specified to meet the dry food 
requirement. Given this harvest requirement, the carbon dioxide, water, and 
nutrient inputs and the oxygen output must be determined for the food producing 
processes. 

The human requirements for oxygen, food, and v/ater have already been given, 
but the corresponding outputs of carbon dioxide and water vapor must be calcu- 
lated. Finally, the waste inputs and characteristics of the vvaste recycler 
must be specified. For details of these calculations see Stahr (1979). 

The recycler is based on an incineration oxidizer proposed by Meissner 
and Modell (1979). The rccycler outputs are chosen to close the material 
balance cycles. This situation corresponds to the ideal case in which all 
waste inputs can be successfully recycled and the CELSS system could survive 
indefinitely as long as each system component performs at its steady state 
operation point without perturbation. 

The next step of the scenario analysis method is the equipment and facil- 
ity specifications. This step along with the survival analysis is essential 
for the break even cost point determination, however, specific details are not 
considered here because they are not needed to illustrate the evaluation process. 

The analysis just completed determines a simplified equilibrium flow 
CELSS model. The steady state mass floi/ values derived in Stahr (1979) are 
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shown in Figure 3. The state equations, each of which expresses mass continu- 
ity for the particular storage involved, are listed in Figure 4, 

A constant growth rate over the plant's 196 day growing cycle is assumed. 
The average biomass holdup (GROW) is half the harvest biomass. So, each day 
of the growing food is harvested. 

To achieve interesting dynamic behavior model development must proceed 
beyond just equilibrium flow considerations. Some examples of the type of 
behavior that should be present in the model include: 

1) The plant growth rate should be affected by environmental conditions. 

2) Required human inputs should be dependent on the level of human 
activity. 

3) Control laws for the recycler and purifier and other flows such as 
the nutrient flow to the growing crop must be formulated. Flows 
needing such laws are indicated with a flow valve indicator in 
Figure 3. 

4) Flov/s from a storage should go to zero as the volume of material 
held in storage approaches zero. 

All these examples and others are included in the model. For instance, 
the plant growth rate is affected by four variables -- the levels of applied 
nutrient and water, the atmosphere carbon dioxide partial pressure, and the 
relative humidity. Some simulations include a fifth variable, a random 
parameter that models a growth rate limit due to some unknown cause. 

A detailed term by term derivation of the model equations is given in 
Stahr (1979). 

4. A GENERALIZED SENSITIVITY ANALYSIS PROCEDURE 

The analysis approach developed by Spear and Hornberger is now presented 
in summary form. (See Hornberger and Spear 1980, Spear 1970a, 1970b for a more 
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detailed exposition). . ,, . , 

4.1 The Mohte-CaHo System Descriptor Analysis Approach 

Consider processes that can be modeled by a set of first order differen- 
tial equations of the form (other models may be handled similarly); 

• ^ ■' dx(.t) = ^(t) = f(2<(.t), e.u.(.t)) 
dt 

where: x(t) = state vector e= uncertain parameters 

u.(t) = inputs, forcing function, time variable functions, 
assumed deterministic here for simplicity. 

Different equations correspond to different CELSS scenarios or different con- 
trol strategies within a scenario. Each element of the parameter vector £ is 
considered to be a random variable with a distribution corresponding to the 
uncertainties in the flows to be expected to or from a CELSS component, or 
to the uncertainty in the required volume of each storage. The parameter 
distributions define an ensemble of models for the scenario. 

For each randomly chosen parameter set e* there corresponds a unique 
state trajectory x(t) and observation vector y*(t). A behavior descriptor 
for the system is defined in terms of y(t). For CELSS this descriptor could 
be survival or nonsurvival of the human occupants for the mission duration. 
More complex multi-possibility behavior categories could be defined. 

Multiple computer simulations are performed for random choices of £ from 
the predefined distributions. Each simulation either results in a behavior, 
B, or not a behavior, B. For each element, of the parameter vector £ 
the multiple runs result in a. set of values ej, associated with B and a set 
associated with F. The subset of parameters that account most strongly for 
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the occurence of B or F must be identified. A large fraction of the parameters 

are often unimportant to the critical behavior. 

A sensitivity ranking can be based on a measure of the separation of the 

cumulative distribution functions F(e|^|B) and F{e|^|F) using the Kolmogorov- 

Smirov two sample test statistic: 

n " sup|S (:x)-S (x)| 
m,n n m ' 

v/here for n behaviors and m nonbehaviors: 

^ = sample distribution corresponding to F(e|^|B). 

= sample distribution corresponding to F(e|^|F). 

Large values of„d^ ^ indicate that the parameter distributions have separated 
and the occurence of B or F is sensitive to the parameter value. In cases 
where induced covariance is small the converse is true for small values of 

^m,n‘ 

A parameter might be strongly coupled with other parameters to affect the 
behavior. Therefore the induced covariance structure must be included in a 
full sensitivity ranking. (See Hornberger and Spear 1980.) 

Upper and lower limit probability bounds for a sample distribution may be 
obtained through use of the statistic 

Dn = sup|S^(x)-F(x)| 

X 

where S^(x) = the sample cumulative distribution. 

F(x) = the true distribution (assumed to be continuous). 

Kolmogorov determined the distribution of D . A confidence band about S (x) 

n n 

that includes F(x) with a given probability can be set up. 

P {D^ = sup|S^(x)-F(x)|> d^} = a 

X 



Given n and a, d may be found in any elementary book on statistics or in Spear 
(1970b). Finally, confidence regions for the true distribution may be found 
from: . / • . ' . . n?; 

P {S^(x)-d^<F(x)<S^(x) + d^for all X} =l-a. ., 

4.2 Model Survival Criteria 

In the model simulations, the behavior, B, is defined as nonsurvival of 
the human occupants. The CELSS must be designed to maximize the probability 
of the survival of the occupants,, F. Failure to survive is evaluated with 
three criteria: 

1) Oxygen. If the human lung alveolar oxygen partial pressure 
ever drops below 35mm Mg the humans perish due to lack of 
oxygen (Billings 1972). 

2) Food. If the humans do not receive their ideal food supply 
a food deficit starts accumulating. Failure occurs if they 
receive the equivalent of no food for a fifty day period. 

Survive if: 

/Tmax 

If a food deficit exists and plenty of food is presently 
available, FM will increase above FMII to start decreasing tlie 
deficit. 

3) Drinking water. Similarly, failure occurs if the humans 
receive the equivalent of no water for a five day period. 



/Tmax 

\ 5l " FHIT^ 

/o 
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Survive if 


Water 


/Tm« 

ii if= U 

/o 


n HMD X .. 
■ HDSS' 


(1 - ^s) dt 


Failure occurs if any of the three criteria indicates failure. 


5. SIMULATIONS 


5.1 Control Lav/s 

Two control strategies are contrasted in the simulations. These strategies 
will be called proportional control and supervisory control. 

The proportional control regulates the flows from storage based solely on 
the amount of matter contained in the storage. The rate laws for three of the 
flows are illustrated in Figure 5. The laws for the other rates, for example 
the waste recycler operation rate (RRTE), are all directly proportional to 
the storage involved (WSTR in this case). 

RRTE = — -g • WSTR • RRAO 

The proportional gains, WSTG in this instance, are all set so that the steady 
state recycler rate is reached when the storage contains 25 days worth of 
material. RRTE is more complicated than the other rates because it contains 
a modulating factor, RRAO, included to model incomplete combustion at low 
atmospheric oxygen concentrations. 

The supervisory control monitors the food and oxygon level and takes 
corrective action if either level starts becoming low. Two steps are taken to 
raise the level 
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1 ) More crops are planted whlcfi produce more oxygen. 

2) The recycler operation rate is decreased Which decreases the 

recycler's oxygen consumptiott. 

V/hen the food storage becomes low the planting ratei is increased. 

. Notice that if the system does become unbalanced, pressure can be applied 
only at certain pbints'such as those above, to attempt a correctibn. However, 
such attempts could conceivably unbalance some other part of the system and 
lead to different failures. 

5.2 Flow Uncertainty Results ■ ' 

. - j vThe . f i rst simulations were conducted with three random parameters relating 
to flow uncertainties. They are: 

1) UHl. Reflects uncertainty in the fraction of the harvest that is 
edible food. 

•2) ULOA. Uncertainty in the human level of activity which affects 
all the human inputs and outputs. 

3) UBOP. A growth rate limit due to some unknown or unmodeled cause 

(i.e., not nutrient, water, carbon dioxide or humidity levels). 

Each parameter was given an independent, rectangular probability density 
function distribution with upper and lov/er limits 25% av/ay from the steady 
state values of 1. Since only upper and lower limits are known, the rectangu- 
lar distribution will provide the least biased specification of the parameters. 

For both control schemes 100 simulation runs of 500 day length were con- 
ducted with different storage volume initial conditions. The results of these 
runs are presented in Table 2. Looking at the results for the Supervisory 
control with 100 day storage it is seen that no failures occurred. Thus, the 
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estimate of the probability of success (survival) is unity. Since 100 simula- 
tions were performed we can say that 

0.96 < probability of survival ^1 

with 95% confidence. (See Clopper'and Pearson, 1934.) 250 days supplies are 
required to make the same statement with the proportional control. ' 

Figure 6 shows state variable plots for one of the 200 day supply propor- 
tional control failures. The humans' perished at 410 days from insufficient 
oxygen. This failure like all the others presented in this paper is a conse- 
quence of the system's dynamics and can not be predicted on the basis of an 
algebraic analysis of flows. 

5.3 Recycler Shutdown 

The response of the model to the unexpected event of a twenty day recycler 
shutdown was studied next. For these simulations the factors UHl, ULOA, and 
UBOP were all set to their steady state value of 1 so that the system could 
operate at steady state in the presence of no distrubances. Storages were 
set to 25 days supplies. 

We might expect system survival since storages exceed shutdown time. 

This expectation, however, is not necessarily correct as can be seen from 
the proportional control results of Figure 7 and the supervisory control 
results of Figure 8. 

With the proportional control, the twenty day shutdown causes a failure 
more than 800 days after the event. The system is never able to re-reach 
steady state. The additional flexibility of the supervisory control allows 
the system to recover without failure. 

The importance of selecting "good" control strategies is illustrated by 
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the improved supervisory control performance in this and the earlier simulations. 

5.4 Storage Size Sensitivity Analysis 

In the previous flow uncertainty simulations we found that 100 days worth 
of supplies were needed to ensure survival. However, some of these storages may 
have been unimportant to survival and been sized much larger than actually neces- 
sary. The sensitivity of individual storage sizes to the recycler shutdown event 
is now studied. The more effective supervisory control is used in the study. 

The storage initial conditions were given rectangular probability density 
functions with the limits indicated in Figure 9. Two hundred fifty simulations 
of 500 day missions were performed which resulted in 103 survivals and 147 fail- 
ures. The cumulative distributions for each storage are shown in Figure 10. 

Table 3' contains the Kolmogorov-Smirov statistics, d^ ^ for each storage. 

The statistic for tv/o storages, NSTR (nutrient storage) and WSTR (waste 
storage), indicate at above the 99% level of significance that the simulative 
distributions fail the null hypothesis of identical distributions. All other 
storages fall below the 99% level. These results indicate, at least in terms 
of a univariate analysis, that the NSTR and WSTR storage sizes are the most 
important determinants of the behavior. 

Note that the distributions for FSTR show that for the initial condition 
ranges used, there are regions such that the survival probability is high for 
cases with small initial food storage. Also, there are regions such that the 
survival probability is low for cases with large initial food storage. How- 
ever, univariate analysis is misleading due to the multiple modes of failure. 

An analysis of the covariance structure should uncover significant couplings 
among the initial condition parameters, if they exist. 
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It may also seem odd that the WSTR initial condition should affect the 
survival probability significantly^ However, the supervisory control gains 
part of its flexibility from being able to decide how fast to recycle the wastes 
in sotrage. This flexibility is lost if only small amounts of waste are present. 

The univariate analysis indicates that the NSTR and WSTR storages are 
adversely affecting the CELSS survival as presently sized. The next step in the 
storage size design process is to increase these storages and repeat the multi- 
ple simulations. - 

NSTR was increased from its previous range of values to a constant value 
of 200 days supply. The large figure was chosen because the NSTR separation 
was so large and the additional nutrient weight is negligible compared to the 
total weight of the CELSS storages. WSTR was increased to the upper limit of 
its previous range. The full initial condition set for this second set of 
multiple runs may be found in Figure 11. 

In a set of 250 simulations absolutely no failures occurred with these 
new initial conditions.. The estimated probability of survival is unity. 

Thus, with 95% confidence we can say that 

P..98 < probability of survival ^1 

for the modified initial conditions when subjected to this particular unexpected 
event . (See Clopper and Pearson, 1934.) With this level of confidence we can 
now conclude that the storage initial conditions as presently sized provide a 
high probability of survival for the recycler shutdown event. 

We might wish to examine a design alternative that does not involve 
increasing the waste storage size. If we want to change any of the initial 
condition probability densities, the entire analysis process must be repeated. 

A multiple run could be done in which all WSTR initial conditions v;ere set 
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equal to zero. New sets of . cumulative distributions would be determined that 
might now indicate different sensitivities caused by the changed system struc- 
ture. : ‘ 

Similarly a different unexpected event could require different storage 
sizing to ensure survival. The total range of possible events, rather than a 
single event, must be defined and considered in a complete CELSS scenario 
evaluation. 

6. DISCUSSION 

An approach to the problem of comparing alternative CELSS of PCELSS sce- 
narios and control strategies within a scenario has been illustrated in this 
paper. We see that once probability distributions have been assigned to all 
uncertain parameters of the system, this method performs a generalized sensitiv 
ity analysis. This enables us to determine the critical parameters or regions 
within the range of the parameters leading to a high probability of survival. 

If the model does not include provision in some form for a possible criti 
cal failure mode of the real system, there is nothing in the sensitivity anal- 
ysis method that will uncover the sensitivity. For instance, microelements 
may build up in the plants or humans or microorganisms might be present that 
could lead to other possible failures. 

The existence of failure modes that are a consequence of the system's 
dynamic interactions and not predictable from static analysis, complicates 
the modeling task in CELSS systems because it is more difficult to determine 
the elements of the system that can be neglected in simulation studies. The 
generalized sensitivity analysis presented here should help in that effort 
as well as find use in control strategy evaluation. 



-15- 

REFERENCES 

Best, C.H., Taylor, N.B., 1963, The Human Body — Its Anatomy and Physiology . 
(Holt, Rinehart & Winston), part 5. p. 296. 

Billings, C.E., 1972, Chapter 2: Atmosphere, in: Bioastronautics Data Book 

NASA.SP-3006, Second edition., 

Clopper, C.J., and Pearson, E.S.,.1934, "The Use of Confidence or Fiducial 
Limits Illustrated in the Case of the Binomial", Biometrika , V. 26, 
pp, 404-413. 

Hornberger, G.M., Spear, R.C., 1980, An Approach to the Preliminary Analysis 
of Environmental Systems. Water Research (in print). 

Incropera, F.P., 1975, Leaf Photosynthesis: The Influence of Environmental 

Variables. J. Environ. Qual., vol. 4, no. 4. 

Lehninger, A.L., 1975, Biochemistry , second edition (Worth), p. 249, 840. 

Lindgren, B.W., 1976, Statistical Theory , third edition (Macmillan), p. 494, 
581. 

Macelroy, R.D., Averner, M.M., 19TT, Space Ecosynthesis: An Approach to the 

Design of Closed Ecosystems for use in Space, NASA Tech. Memo 78491. 

Meissner, H.P,, Modell, M. , 1979, Recycling Plant, Human and Animal Wastes 
to Plant Nutrients in a Closed Ecological System. Am. Soc. Mech. Eng. 
Paper 79-ENAs-29. 

Modell, M. , 1977, Sustaining Life in a Space Colony. Technology Review. 
July/August 1977. 

Modell, M. , 1979, Development of a Prototype Experiment for Treating CELSS 
and PCELSS Wastes to Produce Nutrients for Plant Growth. Research pro- 
posal submitted to NASA. 



-16- . 

Modell, M. , Spurlock, J.M., 1979, Closed-Ecology Life Support Systems (CELSS) 
for. Long-Duration, Manned Missions. Am. Soc. Mech, Eng. Paper 79-ENAs-27. 

Ricci, B. , 1967, Physiological Basis of Human Performance . Philadelphia, Lea 
Febiger, p. 118-9. 

Spear, R.M. , 1970a, Monte Carlo Method of Component Sizing. J. Spacecraft and 
Rockets, Vol. 7» No. 9, Sept. 1970. 

Spear, R.M., 1970b, The Application of Kolmogorov-Renyi Statistics to Problems 
of Parameter Uncertainty in Systems Design. Int. J. Control. Vol.ll, 

No.5. 

Stahr, John, D., "An Approach to the Preliminary Evaluation of Closed-Ecology 
Life Support System (CELSS) Scenarios and Control Strategies", Master's 
Project, U.C. Berkeley, Dec I 1979. 

Tiwari, J.L., Hobbie, J.E., 1976, Random Differential Equations as Models 
of Ecosystems-M. Initial condition and parameter specifications in 
terms of maximum entropy distributions. Math. Biosci. 31:37-53. 

Quinlan, A.V. , 1975, Design and Analysis of Mass Conservative Models of Eco- 
dynamic Systems. Doctoral Dissertation. Dept, of Civil Eng., MIT, 
Cambridge, Ma. 

Quinlan, A.V., Paynter, H.M. , Some Simple Nonlinear Dynamic Models of 

Interacting Element Cycles in Aquatic Ecosystems. Am. Soc. of Mech. Eng. 
Paper 76-Aut-W. 




INPUT 

Gram Person”^ Day” ^ 

Food (Dry) 

600 

Oxygen 

. 900 

Drinking Water 

1 800 

1 

Sanitary Water 

2300 

Domestic Water 

1 6800 


TABLE 1 : Inputs required to support a person in space. 

Adapted from Modell (1977) 



NUMBER OF STEADY 
STATE DAYS STORAGE 
FOR ALL STORAGES 

NUMBER OF FAILl 
(500 day missi( 
"PROPORTIONAL ' 
CONTROL 

JRES IN 100 RUNS 
Dn length) 

SUPERVISORY" 

CONTROL 

50 

39 

22 

100 

21 

1 

0 

150 


0 

200 

6 

- 

250 

. 0 

- 


TABLE 2: Flow uncertainty results 


STORAGE 

*^m,n 

GROW 

.1345 

FSTR 

.1763 

NSTR 

.6561 

AOST 

.1576 

- HSTR 

.1019 

CSTR 

.1567 

WSTR 

i .2510 i 

SDST i 

\ 

— .. 

^ .0965 1 


Acceptance limit for the Kolmogorov-Smirnov tvro sample 
test (Lindgren 1976) 


o< = .01 



103+ IV?' _ 
103 147 " 


.2095 


TABLE 3: 


Kolmogorov-Smirnov two sample test statistic 
for each storage. 103 survivals, 1/f? failures. 
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FIGURE 1 : Simple ecosystem cycle 
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FIGURE 3 I Equilibrium mass flow 
C^SS model. 


HUMANS 


( 9 100 ^ DOMESTIC 
STORAGE 


1-0^--Sjpurifier \-" - - 


Sanitary ^Domestic 19 100 










• ■ ■ ■ -23- •. , ■ . 

^ GROW = BO - 

^^FSTR = .133 ‘ urn • - FM 

^WSTR = .867 - UH2 • + WM - WRSS • RRTE 

^NSTR = NRSS • RRTE - NF 

^SDST = SDM - SOPS • PRTE 

^AOST = OF - OM - OWSS • RRTE 

j^HSTR = HPSS • PRTE + AHF + AHM + AHUS - RRTE - HM - HF 

^ CSTR = CM + CWSS • RRTE - CF 

STATE EQUATION NOMENCLATURE 

AHF Plant atmospheric water output 

AHM Human water output to atmosphere 

AHWS Steady state waste recycler water output 

AOST Atmospheric oxygen storage 

BO Plant net growth rate 

CM Human carbon dioxide output flow 

CSTR Carbon dioxide storage 

CWSS Steady state waste recycler carbon dioxide output 

FM . Actual flow of food to man 

FSTR Food storage 

GROW Growing food crop 

HF Actual plant water flow 

HM Water flow to humans 

HPSS Steady state purifier output 

HSTR Water storage 

NF Actual nutrient flow to crop 

NRSS Steady state recycler nutrient production rate 

NSTR Nutrient storage 

OF Oxygen output rate of plants 

OM Actual oxygen flow to humans 

OWSS Steady state waste recycler output 

PRTE Purifier control rate 

RRTE Recycler operation rate 

SOM Sanitary and domestic water flow to humans 

SOPS Steady state purification rate 

SDST Sanitary and domestic storage 

UHl Edible crop fraction uncertainty 

UH2 Inedible crop fraction uncertainty 

WM Actual human waste production 

WRSS Steady state v/aste recycler input 

WSTR Waste storage 

FIGURE 4: Model state equations and nomenclature 
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FIGURE 5: Flow laws' 

a) Ideal flow until FSTR contains 5 
days supply of food. 

b) Ideal flov/ until HSTR contains 7,4 days supply 
of drinking v/ater for the humans (.5 days 
total water requirement), 

c) Ideal flov/ until HSTR contains 23,6 days supply 
of water for the growing crop (5 days total water 
requirement). 




PROPORTIONAL CONTROL 


Initial Conditions 
(200 days suppli;:s) 


GROW=519027.6 

FSTR= 140880 

NSTR=200 

A0ST=332000 

HSTR=5302040 

CSTR=548000 

WSTR=0 

SDST=0 


Coniments: •UH 1 = 1, 094 

UL0A= 1.224 
UB0P=.7616 

•Failure at 410 days due to lack of oxygen. 


FIGURE 6A: 200 day supply j uncertain flow, proportional control , failure. 








TIME CDAYS> 


TIME CDAYS5 

FIGURE 6C: 200 day supply, uncertain flow, proportional control, failure. 








PROPORTIONAL CONTROL 


Initial Conditions 
(25 days supplies) 

GR0V/=519027.6 
FSTR=17610 
NSTR=25 
A0ST=Zf1500 
HSTR= 662755 
CSTR=68500 
WSTR= 122995 
SDST=Z^77500 


Comments: •Recycler shu-fcoff for 20 days starting at day 10 

•UHljULOA & UBOP all set equal to 1 , 

•Failure at 880 days due to lack of food. 


FIGURE 7A: Recycler shutdown, proportional control, failure. 














SUPERVISORY CONTROL 


Initial Conditions 
(25 days supplies) 

GR0W=519027.6 
FSTR= 17610 
NSTR=25 
A0ST=41500 
HSTR=662755 
CSTR= 68500 
WSTR= 122995 
SDST=if 77500 


Comments: •Recycler shutoff for 20 days starting at day 10, 

•UHl ,UL0Aj&UB0P all set equal to 1, 

•Failure at 880 days due to lack of food. 


FIGURE 8A: Recycler shutdown, supervisory control, survival. Compare 

with figure 7, 



















SUPERVISORY CONTROL 


flNITIAL CONDITION 

LOWER LIMIT 

UPPER LIMIT 

GROW 

.8*519027.6 

1 . 2 . 519027.6 

FSTR- 

.4*17610 

1.2*17610 

NSTR 

.4*25 

1.2*25 

AOST 

1*41500 

1 , 2.41500 

HSTR 

.4*662755 

1.2*662755 

CSTR. 

.4*68500 

1.2*68500 - 

WSTR 

0 

122995 

SDST 

0 1 

^ 477500 


Comments: •I.C,'. s selected from rectangular probability 

density functions. 

•UHljULOA &UBOP all set equal to 1 

•The atmospheric volume is dependent on the AOST I.C. 


FIGURE 9: Stprage size sensitivity analysis initial conditions. 




CUMULATIVE DISTRIBUTION CUMULATIVE DISTRIBUTION 


•40000 


50000 
SROW IC. 


60000 


iOOOO 


15000 
FSTR IC. 


20000 


10 

20 

30 

•40000 

-45000 

50000 


NSTR IC. 


t 

AOST IC. 


FIGURE 10A: 

Cumulative 

distributions. 








CUMULATIVE DISTRIBUTION CUMULATIVE DISTRIBUTION 








SUPERVISORY CONTROL 


p_ 1 

INITIAL CONDITION 

LOWER. LIMIT 

UPPER LIMIT 

GROW 

,8*519027.6 

1.2*519027.6 

FSTR 

,4*17610 

1.2*17610 . 

NSTR 

200 

200 

AOST 

1*41500 • 

1.2*41500 

HSTR 

.4*662755 

1.2*662755 

CSTR 

.4*68500 

1 1.2*68500 

WSTR 

122995 

1 122995 

• SDST 

0 

477500 


Comments: •I.C's selected from rectangular probability 

density functions, NSTR & WSTR increased from 
figure 9. values, 

•UHljULOA & UBOP all set equal to 1 

•The atmospheric volume is dependent on the.AOST I.C 

FIGURE 11 : Increased NSTR & WSTR storage size sensitivity 

‘ analysis initial conditions. 
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Nomenclature 


AHF Plant atmospheric water output 

AHFS Steady state plant respiration fraction 

AHGS Steady state plant lost moisture fraction 

AHM Human water output to atmosphere 

AHMO Human water output to atmosphere due to oxidation 

AHMT Human water output to atmosphere (sv/eating, breathing) 

AHST Water vapor mass per unit volume of atmosphere 

AHWS Steady state waste recycler water output 

ANSI Nitrogen mass per unit volume of atmosphere 

AOST Oxygen storage in total atmosphere , 

BAG Atmospheric carbon dioxide growth rate modulating factor 

BAH Humidity growth rate modulating factor 

BH Available v/ater to ideal water fraction 

BN Available nutrient to ideal nutrient fraction 

BO Plant net growth rate 

BOP Plant total growth rate 

BOPS Steady state total plant growth rate 

BOR Plant respiration rate 

BORS Steady state plant respiration rate 

BOSS Steady state plant net growth rate 

CF Plant carbon dioxide input flow 

CFPS Total plant photosynthesis carbon dioxide flow 

CFRS Plant respiration carbon dioxide flow 

CM Human carbon dioxide output flow 

CMSS : Steady state human carbon dioxide output flow 

C02V. Atmospheric carbon dioxide level modulating factor 

CSTR Carbon dioxide storage 

CWSS Steady state waste recycler carbon dioxide output 

FDEF -Food deficit 

FM Actual flow of food to humans 

FMDF Food deficit flow modulating factor 

FMI . Ideal flow of food to humans (including food deficit) 

FMII Ideal flow of food to humans (not including food deficit) 

FMSS Steady state flow of food to humans 

FRTE Controlled flow level of food to humans 

FSTR Food storage 

FSVL Food survival criteria 

GROS Steady state growing food crop 

GROW Growing food crop 

HDEF Human water deficit 

HUSS Steady state human drinking water flow 

HF Actual plant water flow 

HFGS Steady state plant photosynthesis water flow 

HFI Ideal v/ater flow to plants 

HEMS Steady state plant moisture fraction water flow 

HM Water flow to humans 

HMD Actual human drinking water flow 

HMDF Water deficit modulating factor 



HMDI Ideal human drinking water flow 

HMRT Controlled water rate to humans 

HPSS Steady state purifier output 

HRTA Available water flow to crop 

HSDS Steady state sanitary and domestic water flow 

HSTG Sanitary and domestic control gain 

HSTR Water storage 

HSVL Water survival criteria 

LOA Steady state human level of activity 

LOAA Actual level of activity 

LOAC Carbon dioxide level LOA modulating factor 

LOAF Food deficit LOA modulating factor 

LOAI "Ideal" level of activity 

MFB Supervisory control food level modulating factor 

MOB Supervisory control oxygen level modulating factor 

MOF Supervisory control oxygen level modulating factor 

N Total moles in the atmosphere 

NF Actual nutrient flow to crop 

NFA Available nutrient flow to crop 

NFI Ideal nutrient flow to crop 

NFSS Steady state nutrient flow to crop 

NRSS Steady state recycler nutrient production rate 

NRTE Controlled nutrient flov/ level 

NSTG Nutrient control gain 

NSTR Nutrient storage 

OF Oxygen output rate of plants 

OFPS Steady state total photosynthesis oxygen production rate 

OFRS Steady state plant respiration oxygen consumption 

OM Actual oxygen flow to humans 

OMI Ideal oxygen flow to humans 

OMSS Steady state oxygen flow to humans 

ORTE Oxygen flow to humans modulating factor 

OSVL Oxygen survival criteria 

OWSS Steady state waste recycler output 

P Atmospheric pressure 

PC02 Carbon dioxide atmosphere fraction 

PH20 Water vapor atmosphere fraction 

PN2 Nitrogen atmosphere fraction 

P02 Oxygen atmosphere fraction 

P02L Lung alveolar oxygen partial pressure 

PRTE Purifier control rate 

RH Relative humidity 

RRAO Oxygen level recycler rate modulating factor 

RRTE Recycler operation rate 

SDM Sanitary and domestic water flow to humans 

SOPS Steady state purification rate 

SORT Sanitary and domestic water flow control rate 

SDSG Sanitary and domestic control gain 

SDST Sanitary and domestic waste storage 

SEED Planting rate 

SVL Total survival criteria 

UBOP Grov/th rate uncertainty 
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UHl Edible crop fraction uncertainty 

UH2 Inedible crop fraction uncertainty 

ULOA Random growth rate limiting factor 

VATM Atmosphere volume 

WM Actual human waste production 

WMSS Steady state human waste production 

WRSS Steady state waste recycler input 

WSTG Waste recycler control gain 

WSTR Waste storage 
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